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jy^ . Abstract 

We present a partial proof of van Hoeij-Abramov conjecture about the 
algorithmic possibility of computation of finite sums of rational functions. 
The theoretical results proved in this paper provide an algorithm for com- 
t> ' putation of a large class of sums (1). 
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i> 1 Introduction 
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^ \ The problem of finding sums 

O 



n-l 



S(n) = y £R(k,n) (1) 

k=0 

algorithmically in "closed form" for various terms R(k, n) was intensively studied 
in the last decades (see [1], [2], [4], [10] and references therein). For hypergeo- 
metric R(k,n) which do not depend on n (the so-called "indefinite summation 
problem" ) it was completely solved in [8] ; some further research was done in order 
to make the Gosper's algorithm more efficient. When the terms R(k,n) depend 
on n explicitely (the so-called "definite summation problem") S.A.Abramov and 
H.Q. Le in [2], [4] gave necessary and sufficient conditions for termination of the 
well-known Zeilberger algorithm for rational R(k, n) and general hypergeometric 
R(k, n) thus giving the evidence of incompleteness of the existent algorithmic 
methods in this general case. Recently M. van Hoeij [9], generalizing previous 
observations by S.A.Abramov, proposed a conjecture which (if proved) would 
provide a complete and relatively simple algorithm for finding sums (1) for ra- 
tional R(k,n) = P(k,n)/Q(k,n) e Q[fc, n] if the sum S(n) is again a rational 
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function of n. In this paper we give a partial proof of this conjecture, namely we 
prove that for denominators Q(k,n) which factor in Q(n)[/c] (i.e. over the field of 
algebraic functions of n) as 



and the algebraic functions a s (n) have asymptotics a.i{n) ~ Qn £i with £j < 1 for 
n — > oo, one can algorithmically find out if the sum S(n) is rational and if so find 
it using some simple rearrangements of terms in (1) proposed by S.A. Abramov 
and M. van Hoeij. Such Q(k,n) we will call hereafter "slow denominators". 

The paper is organized as follows. In the next section we explain how one 
can transform the sum (1) in the general case into a finite expression in terms 
of the polygamma functions and algebraic functions. In section 3 we expose the 
transformations of the sum (1) proposed by S.A. Abramov and M. van Hoeij and 
explain their link to the well-known identities for the polygamma functions. In 
that section we also formulate two questions about the possibility of nontrivial 
new identities for the polygamma functions which are crucial for the problem of 
definite rational summation. Sections 4 and 5 give the main theoretical results 
of this paper: the first question is answered completely, for the second we do 
this only for "slow denominators". In the concluding section we discuss some 
algorithmic questions. 

2 Preliminaries 

We always implicitely assume that all terms in (1) are nonsingular: Q(k,n) ^ 
forO<fc<n — l,fceZat least for n large enough. After decomposition of 
R(k,n) = P(k,n)/Q(k,n) into a sum of a polynomial and elementary fractions 
R = P(k,n) + E s (J2ti (fe-afw)O with al S ebraic a «M> ftt(ra) using (2) we can 
represent (1) as a finite sum of a polynomial and of "atomic" sums J2k=o (k-a(n)Y 
with algebraic multipliers (5 st {n). Each atomic sum may be expressed as a differ- 
ence 
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(2) 



s=l 




n — a 



s 



(n)) 



(3) 



^+i)=^)+4 4 



(4) 



where tp( s \z) is the normalized polygamma function 

= (-l) s rf s+1 (logr(^))/(s!^ s+1 ) 
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(cf. [5], [6]). We also use the standard notation ip(z) = d(\ogT(z))/dz = 
T'(z)/T(z) for s = 0. So (1) may be transformed into a finite sum of the form 

£A(n)^( 74 (n)) (5) 

i 

with algebraic fli(n), 7«(n). 



3 Rearrangement of terms and polygamma iden- 
tities 

Obviously one can transform any atomic sum in (3) or directly (1) into an equal 
sum substituting k — > n — k: 

n—l i n— 1 i n—l i 

V - = V - = f-i)T - (6) 

a(n) = n — 1 — a(n) ("Abramov transformation"). This corresponds to the known 
identity 

■0(1 — Z) = l]){z) + 7T COt(7T z) (7) 

and its derivatives for the higher polygammas tp( s \z) after the transformation 
(3). Another transformation proposed by M. van Hoeij splits odd and even terms 
(hereafter we assume without loss of generality t — 1 to avoid cumbersome nota- 
tions): 

2n— 1 y n—l y n—l ^ 

£ (* - «(«)) = £ (2* - + £ (2fc + 1 - a(n)) (8) 

l^ 1 1 l^ 1 1 

~ 2 to (* - « W/2) + 2 to (fc + 1/2 - a(n)/2) ' 

We may choose to split the sum into m parts in a similar way. This transformation 
corresponds to the identity (see [6]) 

1 m ~ 1 / r\ 

i[)(mz) = logm H ip [ z -\ ) (9) 

m ^ V to/ 

and its derivatives. Alternatively one can split the sum into two consecutive 
subsums: 

n—l y my n—m—2 y 

£ (k - a(n)) = £ (k - a(n)) + £ (k + m + 1 - a(n)Y 

This transformation trivializes in the language of the polygamma functions: 
tp{a) — ip(b) = (ip(a) — ip(c)) + (ip(c) —ip(b)). Obviously one can also split off 
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or add a fixed number of terms in the beginning or the end in any sum (1) ex- 
tracting rational functions R(0, n), R(l,n), ... — this corresponds to the basic 
identity (4). 

Now one can formulate the following natural questions: 

Question 1. What finite sums of the form (5) with algebraic (3i(n), 7«(n) are 
in fact equal to rational (or algebraic) functions of n G C (that is for arbitrary 
complex n)? 

Question 2. What finite sums of the form (5) with algebraic A(n), 7i(n) 
are in fact equal to rational (or algebraic) functions of n G N (that is only for 
positive integer n)l 

We show in the next section that the first question may be easily answered: 
the only possible identities of this form valid for arbitrary complex n are those 
obtained from (4) and (9) after substitution of an algebraic function instead of z 
in the arguments, multiplication with another algebraic function and addition of 
several such identities. 

The second question is more difficult. As (7) shows we can find identities with 
extra terms which vanish for integer values of n. In section 4 we prove that there 
is no additional identities in the case of "slow" 7«(n). 

4 Complex Identities 

Lemma 1 Let a(z) be an algebraic function such that: 

a) a{rii) G N for infinitely many rii G N; 
i 

b) lim — > 0; 

c) a(z) ~ f3z for z — > oo and (3 is a positive real number. 
Then a(z) = f3z + 7 for constant (3, 7, (3 G Q, 7 G Z. 



Proof Since lim — > — for some large N G N and we can always assume that 

i->oom No 

Hi is monotonically increasing, there exist infinitely many pairs of integers n^+i, 
Ui s such that Ui s+ i —Ui s < N . Then in turn we can choose among them infinitely 
many pairs with Ui s+ i— Ui s = Ni for some fixed Ni G N. From c) we conclude that 
a'(z) — > f3 for z — > 00 so for the chosen pairs we have \a{rii a+ i) — a{rii a )\ < M . 
Choose infinite subset of such pairs with a(n is+ i) — a(rii s ) = Mi, Mi G N. We 
have two algebraic functions which are equal in infinitely many integer points so 
they are identically equal: a(z) = a(z+Ni)—Mi. Take z G R such that a(z) has 
no singularity on the real axis for z > z so a(z +Ni) = a(zo)+Mi, a(z +2Ni) = 
a(z ) + 2Mi, . . . The algebraic function a(z) is equal to M x {z — z )/N 1 + a(z ) 
in infinitely many poins so they are equal for all z. □ 
Suppose we have a finite sum (1) which is a rational function of n G C. Then 
it has a finite number of poles. The function i]){z) has simple poles at z — 0, 
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z = — 1, z = —2, ... with residues —1; ip^\z) have at the same points poles 
of order t + 1 with coefficients —1 in the Laurent expansions. In order to have 
finite number of poles in (1) almost all must cancel. They can not be suppressed 
by the algebraic coefficients Pi(z) which have only finite number of poles and 
zeros (on an appropriate Riemannian surface). Consider an accumulation point 
of the poles of a given ip^(ai(z)), i.e. a point z e C, such that a^z) — > oo 
for z — > z . Take the highest value of in (5). Then there are several terms 
with the same (t) and the same accumulation point z$. Let the corresponding 
1i P ( z ) ~ c p (z — Zo) £p for z — > z and Eo be the largest exponent among them. Then 
obviously we may consider only the terms with the same value of the exponent 
e in the asymptotics: the smaller values will give density of poles not enough 
for cancellation with this term. Suppose |ci| is the largest absolute value of 
coefficients in this asymptotics (probably we have several such q). Near z make 
the following change of variables: w = a^z), so for z — ■> z , w — > oo. On 
the w-plane in the neighbourhoud of infinity we now consider only the terms 
chosen above, ip® (cti(z)) = ip®(w) has the poles at w = 0, w = —1, w = —2, 
. . . , almost all must cancel with the poles of the finitely many similar terms 
ip^(ai p (z)) = ip^ (Pi P ( w )) j Pi P ( w ) ~ c p w, w — > oo. We easily conclude that 
some of Pi p (w) satisfy the hypotheses in Lemma 1 so they are Pi p (w) = f3 p z + q p , 
Pp ^ Q) Sp ^ Z. All terms in question can be collected in subgroups with 
this property, i.e. the arguments in ip differ only by a transformation ji p (n) = 
Pplio( n ) + Pp £ Q, ^ £ Z.; each subgroup must have separate cancellation 
of almost all poles. Considering now the residues of the poles and using similar 
asymptotic expansions for the algebraic coefficients Pi (z) we conclude that that 
the corresponding coefficients Pi (z) in a subgroup form subsubgroups where they 
differ only by a constant multiplier and ensure separate cancellation of almost all 
poles of this subsubgroup. It may be necessary to split some Pi p {z) into sum of 
other Pi(z). Making another change of variables w — > Mw for some large integer 
M and using (9) we see that for each subsubgroup we obtain a sum of the form 
P(w) J2 P k p ip®(w + n p ) with integer n p ; then using (4), we reduce it to a rational 
function of w multiplied by an algebraic function so it is an algebraic function 
of the initial variable z and may be transformed to it algorithmically. In the 
process of computation one may use the standard Puiseux expansions based on 
the Newton polygon for the algebraic functions. 

Thus we have proved that there are no essentially new identities of the form 
Pi(n)'ip^ (^(n)) = R(n) with rational (and even algebraic) R(n) valid for 
arbitrary complex n, all such identities may be deduced from (4), (9) and for any 
sum (5) we may algorithmically decide if it is a rational (or algebraic) function. 
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5 Integer Identities 



As we have already mentioned the case of identities of the form 

J2fr(n)^ u \l*(n))=R(n) (11) 

i 

valid only for positive integers n is more difficult. First of all we restrict our 
attention to sums of terms Pi(n) (ip^ (n — ai(n)) — ^^(—a^n))) originating 
from (3) and we assume that for n — > oo, ai(n) ~ qn £i with < 1 ("slow 
denominators"). Consider the poles of such terms in the right half-plane Re(n) > 
0. The first term ^^(n — a^n)) obviously has only finitely many of them due to 
our restriction on cnj(n). In order to kill the poles of the second term ip^(— cti(n)) 
let us multiply the sum (1) by Yii sin* i+1 (7r«j(n)). Thus we obtain from the 
original sum (1) a multivalued function without poles or ramification points for 
Re(n) > n . Take a branch of this function for Re(n) > n ; after an integer 
shift we obtain a function f(z) holomorphic in the right half-plane vanishing at 
integers. Using the known asymptotics for the polygamma function in the right 
half-plane we see that even after multiplication by the product of sines its growth 
is bounded as \f(z)\ < Aexp(T\z\ e ) with e < 1. We can use Carlson's theorem 
[7, p. 247]: 

If f{z) is holomorphic for Re(z) > 0, \f(z)\ < Aexp(r\z\) with r < n, and if 
f(n) = for n = 0, n = 1, n = 2, . . . , then f(z) = 0. 

This reduces the situation to that of section 4 so we conclude the absence of 
essentially new identities and algorithmic possibility to decide if a given sum (1) 
with "slow" denominators is rational for sufficiently large integer values of n. 

6 Concluding remarks 

Using the transformations (6) and (8) and their analogs for arbitrary m we can 
make a weaker restriction a.i{n) ~ c(ii e \ £j < 1, but at the moment we see no way 
of treating the general case. Obviously the algorithm that uses our theoretical 
results directly, calculating in the field of algebraic functions Q(n), is far from 
optimal, on the other hand there are known possibilities to avoid some of these 
complications, cf. for example the methods of partial factorization and finding 
dispersion of factors in [1], [3]. 
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